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We show that rain events are analogous to a variety of nonequilibrium relaxation processes in 
Nature such as earthquakes and avalanches. Analysis of high-resolution rain data reveals that power 
laws describe the number of rain events versus size and number of droughts versus duration. In 
addition, the accumulated water column displays scale-less fluctuations. These statistical properties 
are the fingerprints of a self-organized critical process and may serve as a benchmark for models of 
precipitation and atmospheric processes. 



Rainfall and rainfall-related quantities have been 
recorded for centuries 0J^]. All these measurements, 
however, have the disadvantage of low temporal reso- 
lution and low sensitivity. The rain measurements are 
based on the simple idea of collecting rain in a container 
and measuring the amount of water after a certain time. 
The time intervals between readings are typically hours 
or days. Even with the most sophisticated of these con- 
ventional methods, the fine details of rain events can- 
not be captured at all and very light rain might not be 
recorded due to evaporation or insufficient sensitivity of 
the instrument, making it impossible to address questions 
regarding single rain events. 

Recently, high-resolution data have been collected with 
a compact vertically pointing Doppler radar MRR-2, de- 
veloped by METEK |||. The instrument is operated 
by the Max-Planck-Institute for Meteorology, Hamburg, 
Germany at the Baltic coast Zingst (54°43'N 12°67'E) 
under the Precipitation and Evaporation Project (PEP) 
in BALTEX Q. Rain rate, liquid water content, and 
drop size distribution were obtained from the radar 
Doppler spectra, based on a method described by At- 
las ^-0. At vertical incidence, the Doppler shift can 
be identified with the droplet fall velocity. As, in the 
atmosphere, larger drops fall faster than smaller drops, 
spectral bins can be attributed to corresponding drop 
sizes. For a given size, the scattering cross section of 
the droplets can be calculated by Mie theory This 
yields the number density of drops which is proportional 
to the spectral power divided by the corresponding cross 
section. The rain rate q(t) = ^iTiiViVi, where rii is the 
number density of drops of volume Vi falling with veloc- 
ity v%. The detection threshold for rain rates under the 
pertinent operation parameters was q m in = 0.005 mm/h. 
Below this threshold, q{t) = by definition. 

Precipitation profiles up to some thousand meters al- 
titude can be observed. At present, the quantitative re- 
trieval is restricted to rain. Snow and hail can be identi- 
fied from the form of the Doppler spectra but have been 
excluded from the quantitative analysis. The analyzed 
data refer to 250m above sea level and have been col- 
lected from January to July 1999 with 1-min resolution. 



The processes that make a cloud release its water con- 
tent are only very little understood. However, with the 
high temporal resolution of 1 min, single rain events 
can be identified and characterized. Previous work fo- 
cused on the rainfall during a fixed period of time |)]-[ll| . 
What makes the present analysis fundamentally new is 
the identification of a rain event as the basic entity. We 
define an event as a sequence of successive non-zero rain 
rates. Sequences of zero-rain rates in between rain events 
are called drought periods. The event size is defined as 
the released water column in mm, M = ^ t q(t)At, where 
At = 1 min, that is, the time-integral of the rain rate over 
an event. In Fig. 0, the number density of rain events per 
year N{M) versus event size M, is displayed on a double- 
logarithmic plot. In a certain scaling regime, extending 
over at least three decades, the number density of rain 
events obeys a simple power law 

N(M) cx M" 1 - 36 , (1) 



represented by the solid line in Fig. |l| 
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FIG. 1. The number density of rain events per year N(M) 
versus event size M (open circles) on a double logarithmic 
scale. A rain event is defined as a sequence of consecutive 
non-zero rain rates (averaged over 1 min). This implies that 
a rain event terminates when it stops raining for a period of 
at least 1 min. The size M of a rain event is the water column 
(volume per area) released. Over at least three decades, the 



data are consistent with a power law N(M) oc M 
as a solid line. 
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Figure |^ displays the number density of inter- 
occurrence times (drought durations) N(D) between suc- 
cessive rain events. The drought duration is power-law 
distributed 



N(D) oc D~ 



(2) 



implying there is no typical duration of droughts. We 
were not able to detect a lower or an upper cutoff this 
relation. Both the lower end (I min) and the upper end 
(two weeks) still lie within the scaling region of the power 
law. The observed deviation around a period of 1 day is 
related to the daily meteorological cycle. 
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FIG. 2. The number density of droughts per year N(D) 
versus drought duration D (open circles) on a double loga- 
rithmic plot. The drought duration is the time, measured in 
minutes, between two successive rain events. The displayed 
solid line is a power law N(D) oc D~ 142 . The arrow indi- 
cates a time interval of one day. The data deviate from the 
power-law behavior at time intervals corresponding to about 
a day, reflecting the daily meteorological cycle. 

It is compelling, that the distributions of sizes of rain 
events and drought periods are simple power laws. This 
result could prove very useful in relation to drought haz- 
ard assessment or flooding hazard assessment. In order 
to calculate the expected number N(T) of droughts with 
period D > d in a given time period T, we would have 
to integrate N(D). Assuming, for simplicity, that the 
upper cutoff diverges, N(D) = const • D~ 2 , we find 
N(T) = T ■ N(D >d)=T- Sg§|£ • d-° A2 . 

The question of having a reliable water supply is of 
utmost importance. H.E. Hurst |l|,|| posed the following 



problem: How can one design a reservoir so that it never 
overflows or empties? He considered an incoming signal 
q(t) over a time period r. In our case, q(t) is the rain 
rate. The actual level of water in a reservoir (or a river 
) is determined by 



X(t,r) = J>(u)-< g > T )Ai, 



u=l 

where At = 1 min and 



I T 



(3) 



(4) 



denotes the average influx in the considered time period 
r. The water level needed for the reservoir never to empty 
is given by the range 



R(t) — max X(t,r) — min X(t, r) 

Kt<T Kt<T 



(5) 



One can now determine the dimensionless ratio 
R(t)/S(t) as a function of r, where S(t) is the stan- 
dard deviation of the influx q(t) in the period r. For 
uncorrelated random events, this ratio increases as 



R(t)/S(t)cxt h , 



(0) 



where the Hurst exponent H = 1/2. However, Hurst 
GIB discovered that for water level fluctuation in the 
Nile, H w 0.77. Figure || displays the water level X(t, r) 
in a virtual reservoir for the rain data with r = 266, 611 
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FIG. 3. Reservoir level X(t, r) in mm for the entire record 
of duration r = 266, 611 min. The parts of the curve with neg- 
ative slope correspond to dry-periods (droughts) where there 
is no rain, only the mean outflux. The parts of the graph with 
positive slope are periods with rain events. The steepness of 
the line measures the difference between the influxes and the 
average outflux. The range R(t) is indicated with a dashed 
line. 
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Figure ^ demonstrates that R(t)/S(t) oc t h ,H 
0.76 is obeyed over more than four decades with r G 
[10 min, 266,611 min « 6 months]. 
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FIG. 4. Log- Log plot of values of R(t)/S(t) against the 
corresponding period of length r (open circles). The solid 
line is a power law R(t)/S(t) tx t h , H = 0.76. The data are 
consistent with a power law over at least four decades. There 
is a lower cutoff around r ~ 10 min but no upper cutoff is 
apparent. 

It is important to notice that these fluctuations are 
a result of the fluctuating rain rate alone and imply a 
correlation between rain events over the whole temporal 
range studied in this letter. This extends Hurst's result 
H, which he found was valid in the temporal range r 6 
[1 year, 1080 years]. 

It can be illustrated directly that the fluctuations of 
the reservoir are statistically invariant under a transfor- 
mation that changes the time scale by a factor b and the 
level by a factor b H @. In Fig. |, the x-axis of Fig. H 
has been re-scaled with a factor b and the y-axis with a 
factor b H , and the similarity is indeed striking! 
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FIG. 5. Reservoir level X(t, r) in mm for the initial part 
of the record with a duration of r = 74, 743 min. The corre- 
sponding range R(t) is indicated with a dashed line. Com- 
pared with Fig. ^, the x-axis has been rescaled with a factor 
b = 266, 611/74, 743 « 3.57 while the y-axis has been rescaled 
with a factor b H ~ 2.63 ~ 300/115 demonstrating that the 
reservoir level is a self-afhne fractal. 

The power-law number density of rain events is consis- 
tent with a self-organized critical process. The concept 
of self-organized criticality |l5|-|l5[] refers to the tendency 
of non-equilibrium systems driven by slow constant en- 
ergy input to organize themselves into a critical states 
where all scales are relevant. The characteristic feature 
of self-organized critical systems, even if their dynamics 
are incomprehensibly complex, is that the intermediately 
stored energy is eventually released in sudden bursts with 
no typical size. 

A well-known example of such a system is the Earth's 
crust. Currents in the liquid core of the Earth drive 
the crust slowly and fairly constantly. The energy de- 
posited by these currents is intermediately stored in ten- 
sion building up between the tectonic plates and then 
suddenly released in earthquakes. The number of earth- 
quakes per year with a seismic moment 5* exceeding s is 
given by the Gutenberg-Richter law |l(| 

N(S > s) cx s~ B , (7) 

that is, there is no typical size for earthquakes. This 
suggests all the earthquakes have the same physical origin 
and that the Earth's crust is poised in a critical state. 

Avalanches in a pile of grains might also display self- 
organized criticality: when grains are dropped onto a 
pile, one by one, the pile ultimately reaches a stationary 
critical state in which its slope fluctuates about a con- 
stant angle of repose, with each new grain being capable 
of inducing an avalanche on any of the relevant size scales 

0- 

From the perspective of self-organized criticality, rain 
events do not look very different from earthquakes or 
avalanches. If a rain shower, regardless of its duration or 
intensity, is defined as an event, the correspondence to 
avalanches in granular media and avalanches in the crust 
of the earth is striking. The atmosphere is the system 
under investigation and corresponds to the Earth's crust 
or the granular pile. It is driven by a slow and constant 
energy input from the Sun. In particular, water is evap- 
orated from the oceans. The energy is stored in the form 
of water vapor in the atmosphere. It is then suddenly re- 
leased in bursts when the vapor condenses to water drops. 
The power-law distribution of the number density of rain 
events versus size is equivalent to the Gutenberg-Richter 
law for earthquakes and the power-law distribution of 
avalanche size. There is no constant drizzle accounting 
for the constant evaporation but rain events of a wide 
range of sizes. One could imagine having a classification 
of rain events according to their size just as earthquakes 
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are classified according to their position on the Richter 
scale. 

In summary, we found that simple power laws describe 
the number density of occurrence of rain events of a given 
size and drought periods. Moreover, Hurst's analysis 
from the 1950ies on water level fluctuations was extended 
by more than four decades, from a year down to minutes. 
This insight will inevitably inspire new research into the 
modeling of precipitation and atmospheric processes as 
well as serve as a benchmark for existing models and 
might be useful in e.g. drought and flooding hazard as- 
sessment. On a more general level, our analyses show 
that new insights can be obtained from taking the very 
general point of view of complexity and self-organization 
theory. It may serve as an example of how to use this 
approach in situations that seem too complex to be ac- 
cessible to quantitative analysis. 
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